Kondo-lattice screening in a d-wave superconductor 
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We show that local moment screening in a Kondo lattice with d-wave superconducting conduction 
electrons is qualitatively different from the corresponding single Kondo impurity case. Despite the 
conduction-electron pseudogap, Kondo-lattice screening is stable if the gap amplitude obeys A < 
vTkD, in contrast to the single impurity condition A < Tk (where 7k is the Kondo temperature 
for A = and D is the bandwidth) . Our theory explains the heavy electron behavior in the d-wave 
superconductor Nd2- a; Ce :E Cu04. 



I. INTRODUCTION 



The physical properties of heavy-fermion metals are 
commonly attributed to the Kondo effect, which causes 
the hybridization of local 4-/ and 5-/ electrons with itin- 
erant conduction electrons. The Kondo effect for a single 
magnetic ion in a metallic host is well understood 1 -. In 
contrast, the physics of the Kondo lattice, with one mag- 
netic ion per crystallographic unit cell, is among the most 
challenging problems in correlated electron systems. At 
the heart of this problem is the need for a deeper un- 
derstanding of the stability of collective Kondo screen- 
ing. Examples are the stability with respect to com- 
peting ordered states (relevant in the context of quan- 
tum criticality 2 ) or low conduction electron concentra- 
tion (as discussed in the so-called exhaustion problem 3 ). 
In these cases, Kondo screening of the lattice is believed 
to be more fragile in comparison to the single-impurity 
case. In this paper, we analyze the Kondo lattice in a 
host with a <i-wave conduction electron pseudogap 4 . We 
demonstrate that Kondo lattice screening is then sig- 
nificantly more robust than single impurity screening. 
The unexpected stabilization of the state with screened 
moments is a consequence of the coherency of the hy- 
bridized heavy Fermi liquid, i.e. it is a unique lattice ef- 
fect. We believe that our results are of relevance for the 
observed large low temperature heat capacity and sus- 
ceptibility of Nd2- x Ce x Cu04, an electron-doped cuprate 
superconductor—. 

The stability of single-impurity Kondo screening has 
been investigated by modifying the properties of the con- 
duction electrons. Most notably, beginning with the work 
of Withoff and Fradkin (WF)£, the suppression of the 
single-impurity Kondo effect by the presence of <i-wave 
superconducting order has been studied. A variety of an- 
alytic and numeric tools have been used to investigate the 
single impurity Kondo screening in a system with conduc- 
tion electron density of states (DOS) p(w) oc ly T, with 
variable exponent r (see Refs. [6.7.8.9. ld llTlfl2| ). Here, 
r = 1 corresponds to the case of a <i-wave superconduc- 
tor, i.e. is the impurity version of the problem discussed 
in this paper. For r <C I the perturbative renormaliza- 
tion group of the ordinaryi 3 - Kondo problem (r = 0), can 
be generalized^. While the Kondo coupling J is marginal, 
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FIG. 1: The solid line is the critical pairing strength 
A c for T — > [Eq. (|33f) ] separating the Kondo screened 
(shaded) and local moment regimes in the Kondo-lattice 
model Eq. @. Following well-known results^ (see also Ap- 
pendixfS}, the single-impurity Kondo effect is only stable for 
A < L>exp(-2_D/J) ~ T K (dashed). 



a fixed point value J* = r / po emerges for finite but small 
r. Here, po is the DOS for uj — D with bandwidth D. 
Kondo screening only occurs for J* and the transition 
from the unscreened doublet state to a screened singlet 
ground state is characterized by critical fluctuations in 
time. 

Numerical renormalization group (NRG) calculations 
demonstrated the existence of a such an impurity quan- 
tum critical point even if r is not small but also 
revealed that the perturbative renormalization group 
breaks down, failing to correctly describe this critical 
point 9 . For r = 1, Vojta and Fritz demonstrated that 
the universal properties of the critical point can be un- 
derstood using an infinite-?/ Anderson model where the 
level crossing of the doublet and singlet ground states 
is modified by a marginally irrelevant hybridization be- 
tween those state o 10 ^ 1 . NRG calculations further demon- 
strate that the non-universal value for the Kondo cou- 
pling at the critical point is still given by J* ~ f/po, 
even if r is not small^. This result applies to the case of 
broken particle-hole symmetry, relevant for our compari- 
son with the Kondo lattice. In the case of perfect particle 
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hole symmetry it holds that£ J* — > oo for r > 1/2. 

The result J* ~ r/p Q may also be obtained from a large 
N mean field theory 6 , which otherwise fails to properly 
describe the critical behavior of the transition, in partic- 
ular if r is not small. The result for J* as the transition 
between the screened and unscreened states relies on the 
assumption that the DOS behaves as p (w) oc |oj| r all the 
way to the bandwidth. However, in a superconductor 
with nodes we expect that p (to) ~ po is essentially con- 
stant for \u>\ > A, with gap amplitude A, altering the 
predicted location of the transition between the screened 
and unscreened states. To see this, we note that, for en- 
ergies above A, the approximately constant DOS implies 
the RG flow will be governed by the standard metallic 
Kondo result^^ with r = 0, renormalizing the Kondo 
coupling to J = J I (1 — JpohiD/A) with the effective 
bandwidth A (see Ref. 0). Then, we can use the above 
result in the renormalized system, obtaining that Kondo 
screening occurs for Jpo > r which is easily shown to be 
equivalent to the condition A < A* with 

A. = e^ r T K , (1) 

where 

T K =^exp (-JL), (2) 

is the Kondo temperature of the system in the absence 
of pseudogap (which we are using here to clarify the typ- 
ical energy scale for A*). Setting r = 1 to establish the 
implication of Eq. ([1]) for a <i-wave superconductor, we 
see that, due to the ci-wave pseudogap in the density of 
states, the conduction electrons can only screen the im- 
purity moment if their gap amplitude is smaller than a 
critical value of order the corresponding Kondo temper- 
ature Tk for constant density of states. In particular, 
for A large compared to the (often rather small) energy 
scale Tk, the local moment is unscreened, demonstrating 
the sensitivity of the single impurity Kondo effect with 
respect to the low energy behavior of the host. 

Given the complexity of the behavior for a single im- 
purity in a conduction electron host with pseudogap, it 
seems hopeless to study the Kondo lattice. We will show 
below that this must not be the case and that, moreover, 
Kondo screening is stable far beyond the single-impurity 
result Eq. {l}, as illustrated in Fig. Q] (the dashed line 
in this plot is Eq. (D with p = 1/2D). To do this, 
we utilize a the large- N mean field theory of the Kondo 
lattice to demonstrate that the transition between the 
screened and unscreened case is discontinuous. Thus, at 
least within this approach, no critical fluctuations oc- 
cur (in contrast to the single-impurity case discussed 
above). More importantly, our large- N analysis also finds 
that the stability regime of the Kondo screened lattice 
is much larger than that of the single impurity. Thus, 
the screened heavy-electron state is more robust and the 
local-moment phase only emerges if the conduction elec- 
tron ci-wave gap amplitude obeys 

A > A c ~ y / T^D » T K , (3) 
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FIG. 2: (Color online) The solid line is a plot of the Kondo 
temperature 7k (A), above which V — (and Kondo screen- 
ing is destroyed), normalized to its value at A = [Eq. (|14p ]. 
as a function of the d-wave pairing amplitude A, for the 
case of J = 0.3D and p — —0.173. With these parameters, 
7k (0) = 0.0014D, and A c , the point where T K (A) reaches 
zero, is 0.1415 [given by Eq. (133 p ] The dashed line indicates a 
spinodal, along which the term proportional to V 2 in the free 
energy vanishes. At very small A < 2.7 x 10~ 4 D, where the 
transition is continuous, the dashed line coincides with the 
solid line. 



with D the conduction electron bandwidth. Below, we 
shall derive a more detailed expression for A c ; in Eq. ([3]) 
we are simply emphasizing that A c is large compared to 
T K [and, hence, Eq. (J)]. 

In addition, we find that for A < A c , the renormalized 
mass only weakly depends on A, except for the region 
close to A c . We give a detailed explanation for this en- 
hanced stability of Kondo lattice screening, demonstrat- 
ing that it is a direct result of the opening of a hybridiza- 
tion gap in the heavy Fermi liquid state. Since the re- 
sult was obtained using a large- TV mean field theory we 
stress that such an approach is not expected to properly 
describe the detailed nature close to the transition. It 
should, however, give a correct order of magnitude result 
for the location of the transition. 

To understand the resilience of Kondo-lattice screen- 
ing, recall that, in the absence of <i-wave pairing, it is 
well known that the lattice Kondo effect (and concomi- 
tant heavy-fermion behavior) is due a hybridization of 
the conduction band with an /-fcrmion band that rep- 
resents excitations of the lattice of spins. A hybridized 
Fermi liquid emerges from this interaction. We shall see 
that, due to the coherency of the Fermi liquid state, the 
resulting hybridized heavy fermions are only marginally 
affected by the onset of conduction-electron pairing. This 
weak proximity effect, with a small ci-wave gap ampli- 
tude Af ~ ATk/D for the heavy fermions, allows the 
Kondo effect in a lattice system to proceed via /-electron- 
dominated heavy-fermion states that screen the local mo- 
ments, with such screening persisting up to much larger 
values of the d-wave pairing amplitude than implied by 
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the single impurity result^, as depicted in Fig. Q] (which 
applies at low T). A typical finite-T phase diagram is 
shown in Fig. O 

Our theory directly applies to the electron-doped 
cuprate Nd2-sCe x Cu04, possessing both d-wave 
superconductivit y 1 ^ 15 with T c ~ 20K and heavy 
fermion behavior below^ Tk ~ 2 — 3K. The latter 
is exhibited in a large linear heat capacity coefficient 
7 ~ 4J/(mol x K 2 ) together with a large low- frequency 
susceptibility \ with Wilson ratio R ~ 1.6. The lowest 
crystal field state of Nd 3+ is a Kramers doublet, well 
separated from higher crystal field levels^, supporting 
Kondo lattice behavior of the Nd-spins. The supercon- 
ducting Cu-O-states play the role of the conduction 
electrons. Previous theoretical work on Nd2- a; Ce 3 ;Cu04 
discussed the role of conduction electron correlations^. 
Careful investigations show that the single ion Kondo 
temperature slightly increases in systems with elec- 
tronic correlation s 18 ^ 9 , an effect essentially caused by 
the increase in the electronic density of states of the 
conduction electrons. However, the fact that these con- 
duction electrons are gapped has not been considered, 
even though the Kondo temperature is significantly 
smaller than the d-wave gap amplitude A ~ 3.7meV 
(See Ref. We argue that Kondo screening in 

Nd 2 - a; Ce 2 ;Cu04 with Tk *C A can only be understood 
in terms of the mechanism discussed here. 

We add for completeness that an alternative sce- 
nario for the large low temperature heat capacity of 
Nd2- a; Ce a ;Cu04 is based on very low lying spin wave 
excitations 2 -^. While such a scenario cannot account for 
a finite value of C (T) /T as T — > 0, it is consistent 
with the shift in the overall position of the Nd-crystal 
field states upon doping. However, an analysis of the 
spin wave contribution of the Nd-spins shows that for 
realistic parameters C (T) /T vanishes rapidly below the 
Schottky anomaly^ 2 -, in contrast to experiments. Thus 
we believe that the large heat capacity and susceptibility 
of Nd2- a; Ce a ;Cu04 at low temperatures originates from 
Kondo screening of the Nd-spins. 

Despite its relevance for the d-wave superconductor 
Nd2- a; Ce 2 ;Cu04, we stress that our theory does not ap- 
ply to heavy electron d-wave superconductors, such as 
CeCoIn 5 (see Ref. |23j), in which the d-wave gap is not 
a property of the conduction electron host, but a more 
subtle aspect of the heavy electron state itself. The latter 
gives rise to a heat capacity jump at the superconducing 
transition AC (T c ) that is comparable to 7T C , while in 
our theory AC (T c ) < 7 T C holds. 



II. MODEL 

The principal aim of this paper is to study the screen- 
ing of local moments in a e?-wave superconductor. Thus, 
we consider the Kondo lattice Hamiltonian, possessing lo- 
cal spins (Si) coupled to conduction electrons (c kct ) that 



are subject to a pairing interaction: 

H = E ^4a c kQ + g E s < • c Lvi/3 + £W ( 4 ) 

Here, J is the exchange interaction between conduction 
electrons and local spins and £k = £k — ^ with ek the 
conduction-electron energy and [i the chemical potential. 
The pairing term 

[/pair = - E ^kk'Ck"rCL k .l C _k/J.Ck't, ( 5 ) 
k,k' 

is characterized by the attractive interaction between 
conduction electrons £/ kk '. We shall assume the latter 
stabilizes d-wave pairing with a gap A k = A cos 29 with 
9 the angle around the conduction-electron Fermi surface. 

We are particularly interested in the low-temperature 
strong-coupling phase of this model, which can be studied 
by extending the conduction-electron and local-moment 
spin symmetry to SU(N) and focusing on the large- TV 
limit^i. In case of the single Kondo impurity, the large- 
TV approach is not able to reproduce the critical behavior 
at the transition from a screened to an unscreeened state. 
However, it does correctly determine the location of the 
transition, i.e. the non-universal value for the strength of 
the Kondo coupling where the transition from screened 
to unscreened impurity takes place 8 -. Since the location 
of the transition and not the detailed nature of the tran- 
sition is the primary focus of this paper, a mean field 
theory is still useful. 

Although the physical case corresponds to TV = 2, the 
large- TV limit yields a valid description of the heavy Fermi 
liquid Kondo-screened phase 2 -. We thus write the spins 
in terms of auxiliary / fermions as Si • <r a p — > f} a fip — 
<5 Q( g/2, subject to the constraint 

E4/i« = ^A (6) 

a 

To implement the large- N limit, we rescale the ex- 
change coupling via J/2 — > J/N and the conduction- 
electron interaction as C/k,k' — > s _1 C^k,k' [where N = 
(2s + 1)]. The utility of the large- TV limit is that the 
(mean-field) stationary-phase approximation to Ti is be- 
lieved to be exact at large TV. Performing this mean field 
decoupling of H yields 

s 

W = E [& C Ln C km + V (/ km C km + h.c)j + \fl m fkm 
k,m— — S 

s 

- E (A k c_ k _ m c km + /i.e.) +7^0, (7) 

k,m=l/2 

with Eq a constant in the energy that is defined below. 
The pairing gap, A k , and the hybridization between con- 
duction and /-electrons, V, result from the mean field de- 
coupling of the pairing and Kondo interactions, respec- 
tively. The hybridization V (that we took to be real) 
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measures the degree of Kondo screening (and can be di- 
rectly measured experimentally 2 ^) and A is the Lagrange 
multiplier that implements the above constraint, playing 
the role of the /-electron level. The free energy F of this 
single-particle problem can now be calculated, and has 
the form: 



F(V,X,A k ) 



NV 2 NX 



»J]A k A k ^ (8) 



kk' 



N E flfe + A)-^ ktt -Tln(l + e-^- 



k,a=± 



where T = (3 1 is the temperature. The first three terms 
are the explicit expressions for Eq in Eq. ([7]) , and E k ± is 



E k± = 
5 k = 



1 

71 
(A k 



A k 



A 2 + 
A 2 ) 2 



2V 2 + ^ ± 
+ ±V 2 [(& 




describing the bands of our d-wave paired heavy-fermion 
system. 

The phase behavior of this Kondo lattice system for 
given values of T, J and /i is determined by finding points 
at which F is stationary with respect to the variational 
parameters V , A, and Ak. For simplicity, henceforth we 
take Ak as given (and having d-wave symmetry as noted 
above) with the goal of studying the effect of nonzero 
pairing on the formation of the heavy-fermion metal char- 
acterized by V and A that satisfy the stationarity condi- 
tions 



^=0 
dV ' 

dF 

3A=°' 



(10a) 
(10b) 



with the second equation enforcing the constraint, 
Eq. © . We shall furthermore restrict attention to /i < 
(i.e., a less than half- filled conduction band). 

Before we proceed we point out that the magnitude of 
the pairing gap near the unpaired heavy-fermion Fermi 
surface (located at £ = V 2 /X) is remarkably small. Tay- 
lor expanding i?k- near this point, we find 



2 r 



Ek- 



X 
V 2 



(£- V 2 /X- XAl/V 2 



Ak 



1/2 



(11) 



giving a heavy-fermion gap A^k = (X/V) Ak [with am- 
plitude A f = A (X/V) 2 ]. We show below that (X/V) 2 
< 1 such that A /k < A k . In Fig. [31 we plot the 
lower heavy-fermion band for the unpaired case Ak = 
(dashed line) along with ±-E k - for the case of finite Ak 
(solid lines) in the vicinity of the unpaired heavy-fermion 
Fermi surface, showing the small heavy-fermion gap Af k . 
Thus, we find a weak proximity effect in which the heavy- 
fermion quasiparticles, which are predominantly of /- 
character, are only weakly affected by the presence of 
e?-wave pairing in the conduction electron band. 




FIG. 3: The dashed line is the lower heavy-fermion band 
(crossing zero at the heavy-fermion Fermi surface) for the 
unpaired (A = 0) case and the solid lines are ±i?k- for Ak = 
O.ID, showing a small f-electron gap A/k ~ .014D. 




FIG. 4: Plot of the energy bands E+(£) (top curve) and 
£/-(£) (bottom curve), defined in Eq. (113[) . in the heavy 
Fermi liquid state (for A = 0), for the case V = 0.2D and 
A = 0.04D, that has a heavy-fermion Fermi surface near 
£ = D and an experimentally-measurable hybridization gap 26 
(the minimum value of E+ — E-, i.e., the direct gap) equal 
to 2V ~ \?TkD. Note, however, the indirect gap is A ~ Tk. 



III. KONDO LATTICE SCREENING 

A. Normal conduction electrons 

A useful starting point for our analysis is to recall 
the well-known3I unpaired (A = 0) limit of our model. 
By minimizing the correpsonding free energy [simply the 
A = limit of Eq. p|], one obtains, at low temperatures, 
that the Kondo screening of the local moments is repre- 
sented by the nontrivial stationary point of F at V = Vq 



and A = A = V 2 /D, with 



Vo 



! D + H 
2 Po 



■ exp 



2Jp c 



(12) 
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Here we have taken the conduction electron density of 
states to be a constant, po = (2D)^ 1 , with 2D the 
bandwidth. The resulting phase is a metal accommo- 
dating both the conduction and /-electrons with a large 
density of states oc Ao -1 near the Fermi surface at 
ek — A 1 + Vq /Ao, revealing its heavy-fermion character. 
In Fig. [5J we plot the energy bands 



e± (60 = i (; k 



(13) 



of this heavy Fermi liquid in the low-T limit. 

With increasing T, the stationary V and A decrease 
monotonically vanishing at the Kondo temperature 



t k = Hf]y£>2_ Al 2 exp [__L] ; 




(14) 
(15) 



Here, the second line is meant to emphasize that Tk is 
of the same order as the T = value of the /-fermion 
chemical potential Ao, and therefore Tk ^ Vo, i.e., Tk 
is small compared to the zero-temperature hybridization 
energy Vo. 

It is well established that the phase transition-like be- 
havior of V at Tk is in fact a crossover once N is finite-^!. 
Nevertheless, the large- N approach yields the correct or- 
der of magnitude estimate for Tk and provides a very use- 
ful description of the strong coupling heavy-Fermi liquid 
regime, including the emergence of a hybridization gap 
in the energy spectrum. 



B. d-wave paired conduction electrons 

Next, we analyze the theory in the presence of d-wave 
pairing with gap amplitude A. Thus, we imagine contin- 
uously turning on the d-wave pairing amplitude A, and 
study the stability of the Kondo-screened heavy-Fermi 
liquid state characterized by the low-T hybridization Vo, 
Eq. (JT3J) . As we discussed in Sec. HI in the case of a single 
Kondo impurity, it is well known that Kondo screening is 
qualitatively different in the case of d-wave pairing, and 
the single impurity is only screened by the conduction 
electrons if the Kondo coupling exceeds a critical value 



1 



1 



po 1 + lnD/A' 



(16) 



For J < J*, the impurity is unscreened. This result for 
J, can equivalently be expressed in terms of a critical 
pairing strength A*, beyond which Kondo screening is 
destroyed for a given J: 



A* 



D exp [l 



(17) 



[equivalent to Eq. (TTJ for r = 1], which is proportional 
to the Kondo temperature Tk- This result, implying 




FIG. 5: (Color online) Main: Mean-field Kondo parameter V 
as a function of the d-wave pairing amplitude A, for exchange 
coupling J = 0.30D and chemical potential [i = —0.1D, ac- 
cording to the approximate formula Eq. (|31[) (solid line) and 
via a direct minimization of Eq. (JSJl at T = 10~ 4 D (points), 
the latter exhibiting a first-order transition near A = 0.086D. 



that a d-wave superconductor can only screen a local 
spin if the pairing strength is much smaller than Tk, 
can also be derived within the mean-field approach to 
the Kondo problem, as shown in Appendix [5] (see also 
Ref.0). Within this approach, a continuous transition to 
the unscreened phase (where V 2 — > continuously) takes 
place at A ~ A*. 

Thus, calculations for the single impurity case indi- 
cate that Kondo screening is rather sensitive to a d-wave 
pairing gap. The question we wish to address is, how 
does d-wave pairing affect Kondo screening in the lattice 
case? In fact, we will see that the results are quite differ- 
ent in the Kondo lattice case, such that Kondo screening 
persists beyond the point A*. To show this, we have nu- 
merically studied the A-dependence of the saddle point 
of the free energy Eq. (J5J), showing that, at low temper- 
atures, V only vanishes, in a discontinuous manner, at 
much larger values of A, as shown in Fig. [5] (solid dots) 
for the case of J = 0.30£>, fi = -0.1D and T = 1(T 4 D 
(i.e., T/Tk — .069). In Fig. we plot the phase diagram 
as a function of T and A, for the same values of J and 
fj,, with the solid line denoting the line of discontinuous 
transitions. 

The dashed line in Fig.[2]denotes the spinodal T s of the 
free energy F at which the quadratic coefficient of Eq. (JSj> 
crosses zero. The significance of T s is that, if the Kondo- 
to-local moment transition were continuous (as it is for 
A = 0), this would denote phase boundary; the T — > 
limit of this quantity coincides with the single-impurity 
critical pairing Eq. (jTTJ) . An explicit formula for T s can 
be easily obtained by finding the quadratic coefficient of 
Eq. 



1 

J 



E 



tanh£ k /2T s (A) 
2K ' 



(18) 
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with Ek = + ^k' an< ^ w here we set A = [which 
must occur at a continuous transition where V — > 0, as 
can be seen by analyzing Eq. (jlObj) ]. As seen in Fig. [2j 
the spinodal temperature is generally much smaller than 
the true transition temperature; however, for very small 
A — > 0, T S (A) coincides with the actual transition (which 
becomes continuous), as noted in the figure caption. 

Our next task is to understand these results within 
an approximate analytic analysis of Eq. ([8]); before do- 
ing so, we stress again that the discontinuous transition 
from a screened to an unscreened state as function of T 
becomes a rapid crossover for finite N. The large N the- 
ory is, however, expected to correctly determine where 
this crossover takes place. 



conduction electrons. Thus, for T = the particle- 
particle response will be singular. This is the well known 
Cooper instability. For V = we obtain for example 



N D 2 - u 2 
XA (V = 0) = j\o g ^^ : 



(23) 



where we used A as a lower cut off to control the Cooper 
logarithm. Below we will see that, except for extremely 
small values of A, the corresponding Cooper logarithm 
is overshadowed by another logarithmic term that does 
not have its origin in states close to the Fermi surface, 
but rather results from states with typical energy V ~ 

In order to evaluate xa m the heavy Fermi liquid state, 
we start from: 



1. Low-T limit 

According to the numerical data (points) plotted in 
Fig. [5j the hybridization V is smoothly suppressed with 
increasing pairing strength A before undergoing a discon- 
tinuous jump to V = 0. To understand, analytically, the 
A-dependence of V at low-T, we shall analyze the T = 
limit of F, i.e., the ground-state energy E. The essen- 
tial question concerns the stability of the Kondo-screened 
state with respect to a e?-wave pairing gap, characterized 
by the following A-dependent hybridization 



typ 



(19) 



with A typ an energy scale, to be derived, that gives the 
typical value of A for which the heavy-fermion state is 
affected by d— wave pairing. 

To show that Eq. (T9"|) correctly describes the smooth 
suppression of the hybrization with increasing A, and to 
obtain the scale At yp , we now consider the dimensionlcss 
quantity 



XA 



1 d 2 E 



(20) 



that characterizes the change of the ground state en- 
ergy with respect to the pairing gap. Separating the 
amplitude of the gap from its momentum dependence, 
i.e. writing At = A0k, we obtain from the Hellmann- 
Feynman theorem that: 



1 



XA 



2 PQ A \dA/ 

N x-a 

k 



2p A 



km — 



(21) 



For A — > this yields 



XA 



N 
Wo 



du> x 
2^ ^ 

k 



t>lG c 



(k,iw) G cc (-k, -iw) . (22) 



Here, G cc (k,icj) is the conduction electron propagator. 
As expected, xa is the particle-particle correlator of the 



G cc (k,w) 



E- (£„ 



(24) 



where E± is given in Eq. (|13[) and the coherence factors 
of the hybridized Fermi liquid are: 



1 



Inserting G c 
yields 

TV 

XA = -3- 



^/(^-xf + w 2 ) ' 

a -a 



(a - A) + AV 2 



(25) 



(k,cj) into the above expression for xa 



D — u / 4 
■D-u \E+ 



Av 2 u 2 9(EJ) 
E+ + E. 



(26) 

We used that E+ > is always fulfilled, as we consider a 
less than half filled conduction band. 

Considering first the limit A = 0, it holds £L (£) < 
and the last term in the above integral disappears. The 
remaining terms simplify to 



Xa(A = 0) = - 



n r-" „ i 

d£ - 



w 2 ' 



N , D 2 



W 2 



(27) 



Even for A nonzero, this is the dominant contribution to 
Xa in the relevant limit A <C V <C D. To demonstrate 
this we analyze Eq. (|26|) for nonzero A, but assuming 
A -C V as is indeed the case for small A. The calculation 
is lengthy but straightforward. It follows: 



XA 



8 



d r^-^ 2 - 



NX DM 
8 D 8 A 2 



(28) 



The last term is the Cooper logarithm, but now in the 
heavy fermion state. The prefactor X/D ~ T^/D is a 
result of the small weight of the conduction electrons on 
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the Fermi surface (i.e. where £ ~ V 2 /A) as well as the 
reduced velocity close to the heavy electron Fermi sur- 
face. Specifically it holds u 2 (£ ~ V 2 /\) ~ A 2 /U 2 as well 

as£_ (e^yVA)^^(£-^). 

Thus, except for extremely small gap values where 

/ v -D/T K 

A 2 < D 2 i J , xa is dominated by the A = 

result, Eq. ([27| . and the Cooper logarithm plays no role 
in our analysis. The logarithm in Eq. (|27jl is not origi- 
nating from the heavy electron Fermi surface (i.e. it is 

not from £ ~ t- ). Instead, it has its origin in the inte- 
gration over states where < 0. The important term 

4 4 

— m Eq. ([26]) is peaked for £ ~ i.e. where 
J5± (£ ~ 0) = ±V and is large as long as |£| < V. For 
£ ~ holds — — 3W- This P ea ^ at £ ~ has 

its origin in the competition between two effects. Usu- 
ally, u or v are large when E± ~ £. The only regime 
where u or u are still sizable while E± remain small is 
close to the bare conduction electron Fermi surface at 
|£| ~ V (the position of the level repulsion between the 
two hybridizing bands). Thus, the logarithm is caused 
by states that are close to the bare conduction electron 
Fermi surface. Although these states have the strongest 
response to a pairing gap, they don't have much to do 
with the heavy fermion character of the system. It is in- 
teresting that this heavy fermion pairing response is the 
same even in case of a Kondo insulator where A = and 
the Fermi level is in the middle of the hybridization gap. 

The purpose of the preceding analysis was to derive 
an accurate expression for the ground-state energy E at 
small A. Using Eq. (j20|) gives: 



E = E(A = 0) - XAPoA 2 



(29) 



which, using Eq. (|27[) and considering the leading order 
in A <C V and A -C V, safely neglecting the last term 
of Eq. (|28|) according to the argument of the previous 
paragraph, and dropping overall constants, yields 



E 
N 



V A T72 1 A 

-r-^ + V 2 p hx——- 
J 2 D + fi 



PqA 2 ln D 2 



V 2 



(30) 



Using Eq. (flU)) . the stationary value of the hybridization 
(to leading order in A 2 ) is then obtained via minimization 
with respect to V and A. This yields 



V(A) ~ V 



A 2 
WV ' 



(31) 



with the stationary value of A = 2poV 2 , which estab- 
lishes Eq. (fl9|) . A smooth suppression of the Kondo 
hybridization from the A = value Vb [Eq. (|12[) ] oc- 
curs with increasing d-wave pairing amplitude A at low 
T . This result thus implies that the conduction electron 
gap only causes a significant reduction of V and A for 
A ~ A typ cx y/T^D. 

In Fig. [5] we compare V(A) of Eq. (|3T|l (solid line) 
with the numerical result (solid dots). As long as V 



ChoT 




0.0002 0.0003 0.0004 0.0005 



T/D 



FIG. 6: Plot of the low-temperature specific heat coefficient 

% = -WT, for the case of A = 10" 2 D, V = 10~ 1 D, and 
fj, = — 0.1D, for the metallic case (A = 0, dashed line) and 
the case of nonzero d-wave pairing (A = 0.1D, solid line). 
This shows that, even with nonzero A, the specific heat coef- 
ficient will appear to saturate at a large value at low T (thus 
exhibiting signatures of a heavy fermion metal), before van- 
ishing at asymptotically low T < A/ (= A(A/U) 2 = 10" 4 D) 
Each curve is normalized to the T = value for the metallic 
case, 70 a ^ir 2 p V 2 /\ 2 . 



stays finite, the simple relation Eq. (|3Tj) gives an ex- 
cellent description of the heavy electron state. Above 
the small /-electron gap A/, these values of V and A 
yield a large heat capacity coefficient (taking N = 2) 
7 ~ |7r 2 /9 U 2 /A 2 and susceptibility x — 2poU 2 /A 2 , re- 
flecting the heavy-fermion character of this Kondo-latticc 
system even in the presence of a d- wave pairing gap. Ac- 
cording to our theory, this standard heavy-fermion be- 
havior (as observed experimentally^ in Nd2-a;Cea;Cu04) 
will be observed for temperatures that are large com- 
pared to the /-electron gap Af. However, for very small 
T <C Af , the temperature dependence of the heat capac- 
ity changes (due to the d-wave character of the /-fermion 
gap), behaving as C — AT 2 / A with a large prefactor 
A ~ (D/Tk) 2 ■ This leads to a sudden drop in the heat 
capacity coefficient at low T, as depicted in Fig. [5] 

The surprising robustness of the Kondo screening with 
respect to <i-wave pairing is rooted in the weak proximity 
effect of the /-levels and the coherency as caused by the 
formation of the hybridization gap. Generally, a pairing 
gap affects states with energy Ak from the Fermi en- 
ergy. However, low energy states that are within Tk of 
the Fermi energy are predominantly of /-electron charac- 
ter (a fact that follows from our large- N theory but also 
from the much more general Fermi liquid description of 
the Kondo lattice 2 ^) and are protected by the weak prox- 
imity. These states only sense a gap A/k <C Ak and can 
readily participate in local-moment screening. 

Furthermore, the opening of the hybridization gap co- 
herently pushes conduction electrons to energies ~ V 
from the Fermi energy. Only for A ~ V ~ y/T^D will 
the conduction electrons ability to screen the local mo- 
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ments be affected by d-wave pairing. This situation is 
very different from the single impurity Kondo problem 
where conduction electron states come arbitrarily close 
to the Fermi energy. 



2. First-order transition 



the latter we set V = A = in Eq. ([8]), obtaining (recall 



F, 



LM 



N 



AD 



1 A2l 4JT^~- 



-Tln2-T^ln [l + e - ^], 



(36) 



The result Eq. (131|) of the preceding subsection strictly 
applies for A — » 0, although as seen in Fig. in practice 
it agrees quite well with the numerical minimization of 
the free energy until the first-order transition. To under- 
stand the way in which V is destroyed with increasing A, 
we must consider the V —> limit of the free energy. 

We start with the ground-state energy. Expanding E 
[the T — > limit of Eq. to leading order in V and 
zcroth order in A (valid for V — > 0), we find (dropping 
overall constants) 



E 
N 



-4p V 2 In 



A 



16 po 
3 A 



\ 



where we defined the quantity A c 



A c = 4n/d 2 ~ 



/Li 2 exp 



2 Po J 



(32) 



(33) 



at which the minimum value of V in Eq. (|32p vanishes 
continuously, with the formula for V(A) given by 



V(A) 



lAln^ 
2 A 



(34) 



near the transition. According to Eq. (|33[) . the equilib- 
rium hybridization V vanishes (along with the destruc- 
tion of Kondo screening) for pairing amplitude A c ~ 
■\/TkD, of the same order of magnitude as the T = 
hybridization Vq, as expected [and advertised above in 
Eq. ©]. 

Equation (|33p strictly applies only at T = 0, appar- 
ently yielding a continuous transition at which V — > 
for A — > A c . What about T / 0? We find that, for 
small but nonzero T, Eq. (|33[) approximately yields the 
correct location of the transition, but that the nature 
of the transition changes from continuous to first-order. 
Thus, for A near A c , there is a discontinuous jump to 
the local-moment phase that is best obtained numeri- 
cally, as shown above in Figs. and [21 However, we can 
get an approximate analytic understanding of this first- 
order transition by examining the low-T limit. Since ex- 
citations are gapped, at low T the free energy Fk of the 
Kondo-screened (V ^ 0) phase is well-approximated by 
inserting the stationary solution Eq. (j34)) into Eq. ([32]) : 



F K 



1 



A r 



(35) 



for Fk at A —> A c . The discontinuous Kondo-to- local 
moment transition occurs when the Kondo free energy 
Eq. (|35p is equal to the local-moment free energy. For 



where we dropped an overall constant depending on the 
conduction-band interaction. 

The term proportional to T in Eq. (|36[) comes from the 
fact that Fk- = for V = A = 0, and corresponds to 
the entropy of the local moments. At low T, the gapped 
nature of the d-wave quasiparticles implies the last term 
in Eq. (|36[) can be neglected (although the nodal quasi- 
particles give a subdominant power-law contribution) . In 
deriving the Kondo free energy Fk , Eq. (|35p , we dropped 
overall constant terms; re-establishing these to allow a 
comparison to Flm , and setting Flm = Fk, we find 



-p A 2 ln 3 ^ = Tin 2, 
6 P A 



(37) 



that can be solved for temperature to find the transition 
temperature Tk for the first-order Kondo screened-to- 
local moment phase transition: 



T * (A)= 61^2 ln X' 



(38) 



that is valid for A — * A c , providing an accurate ap- 
proximation to the numerically-determined Tk curve in 
Fig. [5] (solid line) in the low temperature regime (i.e., 
near A c = 0.14D in Fig.©. 

Equation (|3"5]) yields the temperature at which, within 
mean-field theory, the screened Kondo lattice is destroyed 
by the presence of nonzero <i-wave pairing; thus, as long 
as A < Tr-(A), heavy- fermion behavior is compatible 
with d-wave pairing in our model. The essential feature 
of this result is that 7k (A) is only marginally reduced 
from the A = Kondo temperature Eq. © , establishing 
the stability of this state. In comparison, according to ex- 
pectations based on a single-impurity analysis, one would 
expect the Kondo temperature to follow the dashed line 
in Fig. H 

Away from this approximate result valid at large TV, 
the RKKY interaction between moments is expected to 
lower the local-moment free energy, altering the predicted 
location of the phase boundary. Then, even for T = 
0, a level crossing between the screened and unscreened 
ground states occurs for a finite V. Still, as long as the 
A = heavy fermion state is robust, it will remain stable 
at low T for A small compared to A c , as summarized in 
Figs. [Hand [2 



IV. CONCLUSIONS 

We have shown that a lattice of Kondo spins coupled to 
an itinerant conduction band experiences robust Kondo 
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screening even in the presence of <i-wave pairing among 
the conduction electrons. The heavy electron state is pro- 
tected by the large hybridization energy ^ > Tr. The 
(i-wave gap in the conduction band induces a relatively 
weak gap at the hcavy-fcrmion Fermi surface, allowing 
Kondo screening and heavy- fermion behavior to persist. 
Our results demonstrate the importance of Kondo-lattice 
coherency, manifested by the hybridization gap, which is 
absent in case of dilute Kondo impurities. As pointed 
out in detail, the origin for the unexpected robustness of 
the screened heavy electron state is the coherency of the 
Fermi liquid state. With the opening of a hybridization 
gap, conduction electron states are pushed to energies 
of order ^JTkD away from the Fermi energy. Whether 
or not these conduction electrons open up a d-wave gap 
is therefore of minor importance for the stability of the 
heavy electron state. 

Our conclusions are based on a large- N mean field the- 
ory. In case of a single impurity, numerical renormaliza- 
tion group calculations demonstrated that such a mean 
field approach fails to reproduce the correct critical be- 
havior where the transition between screened and un- 
screened impurity takes place. However the mean field 
theory yields the correct value for the strength of the 
Kondo coupling at the transition. In our paper we are 
not concerned with the detailed nature in the near vicin- 
ity of the transition. Our focus is solely the location of 
the boundary between the heavy Fermi liquid and un- 
screened local moment phase, and we do expect that a 
mean field theory gives the correct result. One possibility 
to test the results of this paper is a combination of dy- 
namical mean field theory and numerical renormalization 
group for the pseudogap Kondo lattice problem. 

In case where Kondo screening is inefficient and A > 
t/TkD, i.e., the "local moment" phase of Figs. [TJ and 
the ground state of the moments will likely be magneti- 
cally ordered. This can have interesting implications for 
the superconducting state. Examples are reentrance into 
a normal phase (similar to ErRl^B,!, see Ref. [20h or a 
modified vortex lattice in the low temperature magnetic 
phase. In our theory we ignored these effects. This is no 
problem as long as the superconducting gap amplitude 
A is small compared to \JTkD and the Kondo lattice is 
well screened. Thus, the region of stability of the Kondo 
screened state will not be significantly affected by includ- 
ing the magnetic coupling between the /-electrons. Only 
the nature of the transition and, of course, the physics 
of the unscreened state will depend on it. Finally, our 
theory offers an explanation for the heavy fermion state 
in Nd 2 - x Ce I Cu0 4 , where A > T K . 

Acknowledgments — We are grateful for useful discus- 
sions with A. Rosch and M. Vojta. This research was 



supported by the Ames Laboratory, operated for the U.S. 
Department of Energy by Iowa State University under 
Contract No. DE-AC02-07CH11358. DES was also sup- 
ported at the KITP under NSF grant PHY05-5II64. 
APPENDIX A: SINGLE IMPURITY CASE 

For a single Kondo impurity a critical value J* for 
the coupling between conduction electron and impurity 
spin emerges, separating Kondo-screened from local mo- 
ment behavior for a single spin impurity in a (i-wave su- 
perconductor, see Eq. (|I6|) . As discussed in the main 
text, this is equivalent to a critical pairing Eq. (JTTJ) 
above which Kondo screening does not occur. The re- 
sult was obtained in careful numerical renormalization 
group calculations^. In the present section, we demon- 
strate that the same result also follows from a simple 
large- TV mean field approach. It is important to stress 
that this approach fails to describe the detailed critical 
behavior. However, here we are only concerned with the 
approximate value of the non-universal quantity J* . In- 
deed, mean field theory is expected to give a reasonable 
value for the location of the transition. 

Our starting point is the model Hamiltonian 

^ = ekC km C km + + /m/m c km C k'm' 
km m,m', k.k' 

-X)^ kk ' C kT c -ki c -k'i c k'r ( A1 ) 

kk' 

with the corresponding mean-field action S = Sf + Sb + 
Si n t with (introducing the Lagrange multiplier A and hy- 
bridization V as usual, and making the BCS mean-field 
approximation for the conduction fermions): 

S f= l dT H [E4m(»r + ek)c km + f m (d T + A)/ m ] , 
J m k 

S b = J dr^jV^V-XNqo), (A2) 

Sint =J dT J2 (frn^ m V + VWmfm) + £ A k A^ 
mk kk' 

J 

- E ( A kC- k - m c km + c kmC L k _ m A k ), (A3) 

m=l/2 

where the A integral implements the constraint Nqo = 
12 m fmfm, with qo = 1/2. Here, we have taken the large 
N limit, with N = 2 J + 1. 

The mean-field approximation having been made, it is 
now straightforward to trace over the fermionic degrees 
of freedom to yield 



F 



N\V\< 



XNq 



N 



iu — A — Ti(iuj))(iuj + A + Ti(—iuj)) — r 2 (iw)r 2 (ia;) 



(A4) 
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for the free energy contribution due to a single impurity 
in a d-wave superconductor. Here, we dropped an overall 
constant due to the conduction fermions only, as well as 
the quadratic term in Ak (which of course determines the 
equilibrium value of Ak; here, as in the main text, we're 
interested in the impact of a given Ak on the degree of 
Kondo screening), and defined the functions 



i\m = |i/| 2 ]T 



At 



f 2 M = (^) 2 ]T- 



A, 



(iuj) 2 - E 2 ' 



(A5) 
(A6) 
(A7) 



At this point we note that, for a d-wave superconduc- 
tor, r 2 = f 2 = due to the sign change of the d-wave 
order parameter. The self-energy Ti(iu) is nonzero and 
essentially measures the density of states (DOS) Pd{oj) of 
the d-wave superconductor. In fact, one can show that 
the corresponding retarded function Tir(uj) satisfies 



Tir(u>) = |Vf 



dz- 



Pd{z) 



uj + id — z 



(A8) 



with (5 = + , so that the imaginary part T" R (ui) 
—tt\V\ 2 pd(ui) measures the DOS. Writing Tir(u}) 
\V\ 2 G(ui), we have for the free energy 



F 



N\V\ 2 



J 

-N 



XNqo 



(A9) 



f°° dz 

/ — rip (z) tan 

J-oo K 



-\V\ 2 G"(z) 



.z-\-\V\ 2 G'{z))' 
and for the stationarity conditions, Eq. |T 



1 
J 



dz 



n F (z)G"(z){z- A) 



tt (z - A - \V\ 2 G'(z)) 2 + \V\ 4 (G"{z)) 2 ' (A10) 



dz 



n F (z)\V\ 2 G"(z) 



7T {z-\-\V\ 2 G'{z)) 2 + \V\\G"{z)) 



-All) 



which can be evaluated numerically to determine V and 
A as a function of T and A. 

The Kondo temperature Tk is defined by the temper- 
ature at which V 2 — > continuously; at such a point, the 
constraint Eq. (|A11|) requires A — ► 0. Here, we are in- 
terested in finding the pairing A at which Tk — > 0; thus, 
this is obtained by setting T = V = A = in Eq. (|A10|) : 



1 
J 



dz -TTpd(z) 
-D- M T Z 



-Po log 



D + p 



Po, 



(A12) 
(A13) 



where, for simplicity, in the final line we approximated 
p d (z) to be given by 



Pd{w) 



poM/A, for |w| < A, 
po, for \u\ > A, 



(A14) 



that captures the essential features (except for the nar- 
row peak near lu = A) of the true DOS of a d-wave 
superconductor, with po the (assumed constant) DOS of 
the underlying conduction band. 

The solution to Eq. (|A13|) is: 



A, = (D + f i) cxp [1 -1 , (A15) 

Po-J 



showing a destruction of the Kondo effect for A — > A*, as 
V — > continuously, thus separating the Kondo-screened 
(for A < A*) from the local moment (for A > A*) 
phases. 
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